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^ \ Abstract: Flavour changing neutral current decays are a very sensitive test of the 

Q-i! standard model and its extensions. In particular the decay K — > tivu constitutes 

Qh! a clean way to provide constraints, independent of long distance effects. Motivated 

r-| ! by the recent experimental data of the E787 and E865 collaborations and by the 

difference between the standard model prediction and data, we consider in detail new 
^ \ physics scenarios such as the minimal supersymmetric standard model and R-parity 

^ \ violating supersymmetry. We begin with analyzing the impact of new measurements 

on the standard model result obtaining B{K~^ — > n'^uu) = (8.18 ± 1.22) x 10~^^. 
Predictions for other rare kaon decays are discussed, too. Our results allow to improve 
the limits on R-parity violating couplings with respect to previous analyses. 
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1. Introduction 

The rare decay tt+z/z/ is a sensitive probe of quantum effects in the Standard 

Model (SM) and its extensions. We are mainly concerned in this paper with obtain- 
ing limits on the supersymmetric extensions of the Standard Model and in particular 
R-parity violating couplings. R-parity violation allows tree level contributions to 
this decay and can therefore be highly constrained. However the expected smallness 
of the R-parity violating couplings requires that the "usual" supersymmetric loop 
contributions are properly taken into account. Previous studies on upper bounds of 
R-parity violating couplings have been made without taking into account SM and su- 
persymmetric 1-loop contributions. Moreover previous limits do not take advantage 
of the improved experimental results. Another reason to update the previous calcula- 
tions both in the standard and the supersymmetric minimal model is the possibility 
to use the latest data in the flavour sector Q and in related rare decays |]. 
We also take the opportunity to correct a misprint in the neutralino contribution 
contained in the original publications. The impact of this change is quite small: the 
diagrams containing the neutralino contribute to the AF = 1 effective couplings 
together with small down-type Yukawa couplings in the squark mixing in the down 
sector. 

There are numerous classic papers on weak decays and on the analysis of 
the decay tt+z/z/ in the context of the SM and supersymmetry with unbroken 
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R-parity 0. We will follow the improved analyses of 0, ||, |T0| and introduce in 
addition the breaking of R-parity. 

In the first part, we focus on the Standard Model neglecting all new-physics 
effects and we update in this context various related rare kaon branching rates. In 
the second part, we take into account the one- loop supersymmetric contributions 
updating the same rare kaon branching rates. Then in the third part we deduce 
more precise and realistic constraints on R-parity violating couplings involved in the 
Ti^vv decay. Unless otherwise stated all values are at 1-sigma level. 



2. The Standard Model 



Within the standard model the process 
effective Hamiltonian llTll 



TT+z/z/ is governed by the following 



H. 



G 



f 



2a, 



V^vrsin 6^^ 



(2.1) 



where Aj are products of CKM [|lT| matrix elements: Aj = V*gVid. Xt contains the 
top contribution, and X' the charm contribution for flavour /. QCD corrections have 
been calculated at the NLO level 0, [1^ and long-distance effects together with higher 



order electroweak effects are negligible |T^. The top contribution does not depend 
on the lepton flavour since the lepton masses can be neglected with respect to the top 
mass. For the same reason the charm contribution for electrons and muons agrees 
and a function 

Pc{X) = ^^{2Xl + X:) (2.2) 

can be defined, see for the explicit expression. 

The hadronic matrix element for the decay width can be related via isospin to 
the experimentally well known decay 
can be expressed as ||: 



14| and the branching ratio 



BR 



SM 



Im{\t)Xt 
A 



A'^Pc {X) 



i?e(AJ Re{\t 



A 



A 



-X, 



with 



27r2sin^(^,„) 



:2.3) 



[2A) 



0.901 is an isospin violation correction factor |jT3 and A = IV^^I in the Wolfen- 



stein parameterization of the CKM matrix |15]. The error in neglecting effects of 
higher order in A in the improved Wolfenstein parameterization is of the order of 
0.8%. Therefore we can safely neglect these corrections without any significant loss 
in the predictions. 
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The branching ratio for the decay —>■ Tt^e^Ve has recently been measured 
with high statistics by the E865 collaboration [^]. Their result, 

BR{K+ 7r°e+z/e) = (5.13 ± 0.15) x 10"^ 

differs considerably from the most recent value of the Particle Data Group [|16| , 

BR{K+ 7r°e+z/e) = (4.87 ± 0.06) x 10'^ 

which does not include yet the above mentioned result. We will use for our analysis 
an average value for this branching ratio, where we take into account the Particle 
Data Group fit as well as the E865 result: 

BR{K+ 7r°e+z/e) = (5.08 ± 0.13) x 10"^ 

This make the central value increase by 4.4% and gives a larger error on all branching 
rates. Then, k+ = (7.97 ± 0.20) 10"^ . 

The CKM parameters are taken and updated from the recent fit of Stocchi |jl| 

\Vcb\ = 0.04135 ± 0.0007 
A = 0.22385 ± 0.00355 
p = 0.190 ±0.046 
f] = 0.3485 ±0.0275 

These values have been obtained with the most recent value of the top quark mass, 
i^tif^t) = 168.1 ± 4.1 GeV in the MS scheme, corresponding to the experimental 
value for the pole mass of = 178.0 ± 4.3 GeV 



For the charm and top contribution we agree with Ref. |I8| and we use respec- 
tively : 

Pc{X) = 0.39 ± 0.07 
Xt = 1.529 ± 0.042 . 

All these values can be combined to give our standard model prediction for the 
branching ratio of ^ n^uu : 

BR^^^ = (8.18 ± 1.22) X 10"^^ (2.5) 



The central value differs slightly from the recent prediction by Buras et al. |T8| due to 



differences in the CKM parameters resulting from different fits. The fit by Stocchi [|I| 
we used includes some new E865 data and has been updated using the MS top mass 
given above. The error on the branching ratio is around 15%, mostly due to the 



For the relation between the pole mass and top quark mass in the MS scheme see, e.g., Ref. [17 
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charm sector and to p. Buras et al. |18] pointed out that the error in the charm 



sector resulting from the renormahzation dependence of the function Pc{X) could 
be reduced to 0.03 within a NNLO level calculation of QCD corrections. Then the 
main uncertainty in the charm sector arises from m^. 

The theoretical prediction is still compatible with the recent experimental result 
for K+ TT+z/z/, BR = (1.47 tH) 10"^° 0, even if the predicted central value 
is half the observed value. It is too early to speculate about the presence of new 
physics contributions as the experimental error is still large, but it is interesting to 
note that possible new physics effects should be of the same order as the SM ones in 
order to get the measured central value. From these considerations one can conclude 
that the rare decay — >■ 7t~^i>u is likely to play a major role in the future both for 
constraining and discovering effects beyond the standard model. 

The effective Hamiltonian, Eq. |2.1| , governs other related rare kaon decays [|l^ 



too. We can thus easily obtain predictions for the branching ratios of these decays: 
BR''\Kl ^ vrVz/) = KL ^MM^^ ' = (2.90 ± 0.54) x lO"" (2.6) 



BR^^{Kl ^ vrVz/) = Rs 
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= (0.57 ±0.09) X 10"^^ (2.7) 
BR^^^Kl 7r°e+e-) = (3.7 ± 1.2) x 10"^^ (2.8) 



where 



= = (3.49 ± 0.09) X 10"^ (2.9) 

T{K+) rK+ 

Kg = rJ^^^^ = (6.04 ± 0.16) X 10"^ (2.10) 



and = 0.944 is the isospin correction for WA 



Note that the process K'f^ —>■ 7i^e~^e~ is not determined entirely by the short- 
distance contributions governed by the effective Hamiltonian, but there are important 
long-distance effects. These have been included in our estimate for the branching ra- 
tio following the recent work in Ref. [|1^, Using the average value for BR{Ke3^), 
the result does not change significantly with respect to that obtained in Ref. p]|]. 
In contrast to the decay — > n'^uu, experimentally no process has been observed 
so far corresponding to these decays. Hence on the experimental side there exist 
only upper bounds on these branching ratios which are orders of magnitude above 



the present theoretical predictions These decay processes are therefore not well 
adapted to further constrain new physics scenarios. We only mention these branch- 
ing ratios in view of future data as theoretical predictions can be easily obtained on 
the same footing as for n^uu. 
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3. Supersymmetry 



In supersymmetric extensions of the standard model a large number of particles 
carries flavour quantum numbers and can therefore contribute to flavour changing 



neutral current processes [22|. The AF = 1 decay n^uv can be used in a 

model-independent way to constrain new physics effects There are no SUSY 
contributions at tree level, they start only at one-loop order. At first sight the super- 
symmetric contributions thus seem to be small, but as pointed out in Ref. they 
can be of the same order of magnitude as the standard model ones. At one-loop order 
supersymmetric contributions to FCNC processes can be grouped into three classes: 
the Higgs and W/quark exchanges, which contain as a subclass the SM contributions; 
chargino and neutralino/squark exchanges; and gluino/squark exchanges. The dom- 
inant supersymmetric contributions are given by chargino/squark diagrams. In the 
following part of this section we will detail these contributions following Refs. |l^ 



We will throughout this work use the nomenclature "SUSY contributions" for all 
supersymmetric contributions without any R-parity violating ones as the latter will 
be discussed separately. 

The determination of the supersymmetric contribution is based on the same 
principle as in the SM case. There are two classes of diagrams: penguins giving 
rise to an effective Z sd-veiiex and box diagrams. The effective Hamiltonian can be 
written in the same way as in Eq. |2?l| with Xt replaced by Xnew = rKe~''^'^Xt 
tk and 6k measure the amount of new physics and are functions of masses of new 
particles and new couplings. New physics effects proportional to Ac are included in 
the new form of Xt. The branching ratio for the decay n^uu can then be 

expressed as 

= S^ir^Xt I p RtY + (A^P,(1 - 

+2 I Kb |2 (1 - ^yX^P.XtrKRtCosiP + eK)'^ (3.1) 

where p = I — Rt cos (3 and fj = Rt sin [3. 

In writing the effective Hamiltonian in that way, we made some assump- 
tions. First of all we assumed that only one dimension six operator contributes: 
SLl^di i^i7/i^i- At this stage we also assume unbroken R-parity. The field content 
will be chosen minimal (MSSM like). Xnew can then be written in the following form: 

Xnew = Xj'^^^'* + Xh± + X(j + Xfj , (3.2) 

showing the separate contributions in the loops of the standard model particles, 
the charged higgses and up-quarks, charginos and up-squarks (and charged sleptons 
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in box diagrams), neutralinos and down-squarks (and sneutrinos in box diagrams), 
respectively. The gluino contribution (boxes and penguins) has been neglected. In 
the appendix we recall all the formulae needed to compute the SUSY contribution to 
the d svv transition. We have corrected the misprints in X^j and Xh± in Ref. 



originally noticed by the authors of Ref. |T0|. Note also a misprint in the neutralino 



penguins in the expressions in Ref. which has been corrected here (cf. appendix). 

Due to the complicated flavour structure of supersymmetric theories, it is often 
useful to diagonalize the mass matrices perturbatively via the mass- insertion approxi- 
mation |23 . In most cases it is sufficient to perform the single mass- insertion approxi- 



mation. Colangelo and Isidori [TU| have shown, however, that for the chargino/squark 
diagrams, it is necessary to go beyond the single mass- insert ion approximation in or- 
der to take into account all possible large effects. Introducing the R-parameters 
specifying the different mass insertions. 







dL 














V - ^^'^ 



B^st =i^^^^^Vu='^^Vuni.,, (3.5) 
^ « XtmtmuL AtTTLt 



the chargino/squark and the neutralino/squark contributions to Xnew can be written 
as [PlUOl: 



X^ = NRl,^ (3.8) 



rn?- {rn\^) are thereby the squared masses of the down (up)-type squarks of the first 
two generations and (^i))ij(("^^)jj) denotes the corresponding off-diagonal elements 
of the squark mass matrices. The quantities Ci and are listed explicitly in the 
appendix. 

The detailed flavour structure of SUSY models is unknown and the R-parameters 
are thus complex, a priori unknown, quantities. It is however possible to derive upper- 
limits using various experimental resultsf^, |TU|, Since the analyses of [^, |T^] the 
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experimental situation has not changed significantly. Therefore our limits^ 

i?f . -< (-105 - 562) 



's,d, - I ""'^500 GeV 



RY... I < 41 



' ^500 GeV^ 

are very close to their limits. The main changes are due to the new fits to the CKM 
matrix elements (central values). Considering the errors in the CKM matrix elements 
and in the top quark mass, these limits become slightly larger. 

Unfortunately the results for the branching ratio is very sensitive to the SUSY 
parameters (masses and couplings) and only the future data of the new generation 
experiments will give a significant improvement on these quantities. However, it is 
possible to estimate the order of magnitude of the SUSY contributions. The authors 
of Ref. P found for rx and 9k the typical ranges 

0.5 <rK < 1.3, -25° < 9k < 25° . (3.9) 

by varying all SUSY parameters within the bounds allowed by experimental con- 
straints. This analysis is not restricted to some specific model but it is valid for a 
general supersymmetric extension of the standard model provided that new physics 
contributions to the tree level decay —>■ Ti^e^Vf. and contributions due to other 
operators than s'll'^di i^^ltii^L can be neglected. One remark of caution is in order 
here: these ranges are not exclusive but they only indicate the most probable values. 
Our updated analysis agrees with this statement slightly enhancing the probability 
for the values lying within the above range. 

Varying tk and 9k within these ranges makes at most a change of ~50% for the 
branching ratio of K'^ tt^uu. We obtain 

BRS^'y{K+ Ti+uu) = (8.18 ttil) X 10"^^ • 

The central value thereby corresponds to the SM value, i.e. = 1 and 9k = 0. 
Allowing for wider, less probable ranges, we can set an upper limit on the central value 
of approximately 1.55 x 10~^° (with tk = 1.5 and 9k = —30°). Hence it is possible 
to saturate the experimental central value with supersymmetric contributions alone 
without any R-parity violating couplings. 

In the same way as in the standard model case we can deduce branching ratios 
for other rare kaon decays: 

BR^^^y^Kl n%u) = (2.90 ^l^g) x 10"^^ (3.10) 

BR^^^^^^K^-^ n%iy) = (0.57 toU) x 10"^' (3.11) 

^for the first two limits, thcsign means that the upper hmit of the real part (resp. the 
imaginary part) of the R-parameters is the real part (resp. the imaginary part ) of the r.h.s. 
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For the branching of TT^e^e in new physics scenarios see |20|. The range 

given corresponds thereby to the range, Eq. pl9| , for the SUSY contributions and the 



central value corresponds to the SM value. We would like to emphasize that these 
values can only be an order of magnitude estimate since the precise value of the SUSY 
corrections is not known. The range, Eq. p.9| , only gives the most probable value, but 
is in no way exclusive. In addition, we did not perform a really systematic analysis, 
including for instance the errors on CKM parameters. However, these results show 
that the SUSY corrections can be of the same order as the SM ones and that they 
should thus be taken into account for an analysis of the constraints imposed on R- 
parity violating couplings. The latter will be discussed in more detail in the next 
section. 



4. The R-parity violating contribution 

This section is devoted to a discussion of the contributions to the process 

TT^vv arising from R-parity violating couplings. Allowing for R-parity violation in a 

supersymmetric extension of the standard model a superpotential of the form 

W = \,kUL,E, + K^,UQ,Dl + A;;,[/r^P^ (4.1) 

has to be considered. The first two terms of this superpotential violate lepton number 
and the last one violates baryon number. There are rather stringent limits on the 
simultaneous presence of lepton and baryon number violating couplings since they can 
mediate rapid proton decay. In this paper we will concentrate on the A^^^ couplings 
since they induce tree level contributions via squark exchanges to the decay — > 
TT+z/z/ (cf the diagrams shown in Fig. |lD. 

We note that there is an ambiguity in defining the couplings because of an 
ambiguity in defining the lepton and Higgs fields (for a more detailed discussion of 
that point see Refs. p5|, ^). We will take the basis as defined by the superpotential 
{Super- CKM basis), Eq. [4.1| , and assign the CKM matrix to the up- type squarks. 

Including the R-parity violating processes, the branching ratio for the rare decay 

— s> n^uu can be written in the following form: 



R r?( ^ ^+i,r,\ — ^+ I \^ I ^modcl J ^ |2 , \^ I ^bl 

BK{K yv)- — y2_^ I + 4^ ^200 GeV)^ ' + ^ 16^1200 GeV)^ 

(4.2) 

Throughout our analysis we assume (7™°'^^^ to be either Cf*^ = XcX^. + XtXt or 
QSUSY ^ y^^x[ + XtXnen^. The factor "fc" is given by 

k = g^"^"^^^ = 8.88 X 10-« GeV-^ (4.3) 
V^7rsin2(^^) 



|2 
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A' 



2k 



X'n 




Figure 1: R-parity violating diagrams contributing to the process K'^ 



Errors are below 1% and are therefore negligible compared to other errors. The 
R-parity violating couplings are contained in the e^j which are defined as 



E 



X * X 

2 

m - 

driji 



^^^i^ j (200 GeV)' 



2 

m- 



(4.4) 



Note that we give our constraints for degenerate squarks masses of 200 GeV as the 
recent lower bounds on squarks have increased ||T6[. X'^jj^ are complex parameters 
such that the phase of e^j is a priori not known. 

In contrast to the standard model and the SUSY contributions, R-parity violat- 
ing couplings can induce processes with a neutrino and an antineutrino of different 
flavour in the flnal state. This leads to the last term in Eq. [4.2| . Obviously no 
interferences occur with the standard model/SUSY contribution in that case. The 
R-parity violating processes with the same neutrino flavour in the flnal state, how- 
ever, interfere with the SM/SUSY contribution as can be seen from the flrst term in 
Eq. [4.2|. The resulting contribution to the branching ratio is: 



BR 



Int 



12A2 A;(200 GeV)- 



model 



(-11, 



(4.5) 



5. Full analysis and constraints on the A' 



Within this section we will discuss different bounds on the arising from the 
process K'^ ■n^pu. Of course it is not possible to establish bounds on each 
of the 27 couplings separately. Under some assumption it is possible to constrain 
certain combinations of couplings. We will perform our analysis in three steps. 
First we will neglect all contributions except the tree level R-parity violating ones 
in order to estimate their order of magnitude in comparison with the SM and SUSY 
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contributions. The "pure" R-parity violating branching ratio is given by 

^^"^ ^ 48 A2 fc2(200 GeV)4 5Z I I' • (^-l) 



Comparing this with the experimental value it is therefore straightforward to derive 
an upper bound for the sum of all 6,^ : 

I eij p< 5.6 X 10"^° . (5.2) 

With squarks at 100 GeV, we have ^j^- | e^j p< 0.35 x 10^^°. This bound is much 
lower than the upper bounds 2.3 x 10~^ and 2.3 x 10~^° obtained respectively by 
Choudhoury and Roy a few years ago and recently by the authors of Ref. [p^] , 
both using degenerate squarks masses of 100 GeV. On the one hand this shows how 
much data on n^uu has improved recently, but on the other hand this implies 

that the one-loop SM or SUSY contributions should be properly taken into account 
in order to obtain a realistic limit as they are of the same order as the possible 
R-parity violating tree-level ones. 

From our previous discussion, we have drawn the conclusion that SUSY has to 
be included in the analysis of n^vv. In the next step we will therefore include 

the SUSY contribution. Note that this contains the standard model contribution as a 
special case {tk = 1 and 6k = 0). We will therefore not discuss the latter separately. 
An analysis of the resulting bounds including only the SM contribution has recently 
been performed for some special cases in Ref. [ |30|| . 

Since we aim to obtain an upper-bound on the R-parity violating couplings, we 
assume the SUSY contribution (which already includes SM) to be minimal (cor- 
responding to Tk = 0.5 and 6k = 25°) in order to allow for the largest possible 
contribution from R-parity violating terms. To simplify the discussion we will, for 
the time being, neglect interferences (cf. section ^). The upper bound for the sum 
of the eij can then slightly be improved compared to the case without any SUSY 
contribution given in Eq. p.2[ We obtain 

^ I eij p< 4.45 X 10^^° . (5.3) 

This limit can be translated to an upper bound on the product of two couphngs by 
naively setting all the couplings to zero except one product. This procedure results 
in the following bounds: 



2 

m- 



2 

m-. 



< 2.11 X 10"^ (5.4) 

< 2.11 X 10"^ (5.5) 
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Cti 




Cr 


P 




SM 


d (2.03 X 10-5) 


1 


0.89 


5.79 X 10-6 


2.48 X 10-1° + af /4 + [3^4 




q(1.17x 10-5) 


1 


0.81 


5.49 X 10-*^ 


4.47 X 10-10 + J2^ «*V4 + /?V4 



Table 1: Values of the constants appearing in Eq. 5.6 



Of course, more realistic and precise constraints should take into account inter- 
ferences. This, however, makes the extraction of upper bounds harder and no simple 



bounds in the sense of Eqs. |5.3| , ^.4| , |5.5| can be given. In the following part we will 
assume that only final states with the same neutrino flavour occur. Thus only eij 
with i = j has to be taken into account. The general equation verified by the en can 
be written in the following way: 

E (^^M + y)' + E (Me..) + f)' = i?% (5.6) 



Taking only one of the en nonzero, this equation describes a circle in the complex 
plane, whose parameters are listed in Table |l| for the case of the standard model {rx = 
1 and 6k = 0) and for the "minimal" SUSY contribution, corresponding to rx = 0.5 

48\2 1.2 

and Ok = 25°. The radius R contains the term — — —{BW'p - BR^^) and 

/€+(200 GeV)-4^ ' 

the shifts a and (3. The shifts depends on the CKM inputs and on the loop-functions 
Xt and Xl The the corrections of a due to the implicit dependence on the 

lepton masses mj. 

The resulting constraints in the complex plane on en are displayed in Fig. ^ 
Constraints on €22 and €33 can be obtained in the same way and are of the same 
order of magnitude. To have an idea of the influence of the interferences, we may 
choose the point of coordinates (Re(eii)=-2, Im(eii)=-2) on the SUSY circle of Fig. 2. 
It is approximately the point which gives the maximum value for | en |. We have 
I en 1= 2.8 x 10-^ and then we deduce, in the same manner, limits on products of 
RPV couplings'^: 



< 2.8 X 10-5 (5.7) 

A * A' 

inl in2 \ ^ o o 1 n— 5 



2 



2 

m-; 

driL 



< 2.8 X 10-^ (5. 



These are 30% bigger than the ones in Eqs. |5.4|, [5.5| and so this numerical example 



shows that interferences may have a significant influence. 



^the case \=t gives slightly lower limits due to the correction factor Ct — 0.81, but the limits 



5.7, 5.8 can be used for the 3 flavours. 



- 11 - 



-3 



-2 



-1 
Re 



1 



Figure 2: Allowed region for i?e(eii) and Im(eii) in units of 10~^. In order to compare 
with the very recent constraints on these R-parity violating couplings [^], remember that 
we take 200 GeV as reference value for the mass of the squarks. 



6. Conclusions 



We have investigated the decay n^uu and related rare decays as a probe 

of physics beyond the standard model. As a starting point we have obtained 
the standard model value using updated experimental values for all the relevant 
parameters. We found a slightly bigger branching than other recent estimates: 
B = (8.18 ± 1.22) X 10^^^. Furthermore we have analyzed the supersymmetric 
contributions in the mass-insertion approximation and corrected a minor misprint in 
the neutralino contribution that was present in the existing literature. The main con- 
cern of this paper was to obtain stringent limits on the R-parity violating couplings. 
Assuming that the process — >■ tt'^uu is governed entirely by RPV processes, 
the bounds on RPV couplings can be lowered with respect to previous analysis (cf. 



Refs. p9| , |26[| ) due to recent data. Recent experimental limits, however, indicate that 
SM and SUSY contributions, which can give up to 50% of the SM one, should be 
taken into account in order to establish limits on the RPV couplings. We performed 
an analysis including all these effects and we have shown that interferences can have 
an effect of the order of 30 % on the bounds on RPV couplings. 
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on R-parity violation. Feynman diagrams are drawn using Jaxodraw |31 



A. Explicit expressions of the supersymmetric contributions 

We only list the final formulae, more details and explanations (on the basis of the 
fermions and sfermions fields, Feynman graphs) can be found in Refs. 0, |10|. The 



misprints of ^ have been corrected and functions are written in the notation of [|T^ 
We based our calculations on |5B |. 



General notations 



Xij denote ratios of squared masses (for example : 



• j(xi,..,x„) and k(xi,..,x„) are the loop functions defined in (singularities 
of these functions in the case of equal arguments become derivatives). 

Charged Higgses contribution 

Chargino contributions 

• Qt = IS the top-quark couplmg, 

\/2Mw sin P 

• U and V are 2x2 unitary matrices that diagonalize the chargino mass matrix: 

= diag{Mc^,Mc^) V, M^^ > (A.2) 



They can be found in explicit form in [Rl] [RH] 
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1-RR contribution 



c'rr = E i^^')^r + ) (A-4) 

= I ^m- Jfe., ^.f., ^e.fj 4^^2C/a4 (A.6) 



2 

2-LL contribution 



i;i=i,2 



(A.8) 

3-LR and RL contribution 



"2 j [xiu]^ , Xju^ , ) L'ji ?7|j 



2mii^ 



(^«^)r = -|;f^4^^2X (A.14) 
(Crl)j^''' = -^^^ j(a;is^, a;^-^^, a;t^ni, a;e^«L)^ij^i2^/i2^/2j (A.15) 
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Neutralino contribution 



The neutralino contribution is potentially large since -Rf^rf^ can be large such that 
it should be taken into account. Compared to Ref. 0, we corrected a minus sign 
in front of the function k in expression for the penguin diagram. An overall missing 
factor (-1/2) has been corrected too. The first two diagrams in figure 7 in Ref. 
involving a squark-squark-Z vertex are cancelled by self-energy corrections. The two 
other contributions give: 

^= E ((^)nr+ (A^)r) (A.16) 

n;m=l.A 

-^k(x„,-^,x^,-^, 1) w^{d,)[w:,wl, - w:,wl^]w:;,{d,) (a.i?) 

(A.18) 

W is the unitary matrix which diagonalizes the Neutralino mass matrix : 

= W^.dtag{Mj^.M^. M^.M^YW , M^^ > (A.19) 



Mm 



/ Ml 

M2 
—Mz sin cos [3 Mz cos Qw cos (3 
y Mz sin Qw sin jd —Mz cos Qw sin (3 

Y{x) 



—Mz sin Qy/ cos (3 Mz sin Qy/ sin j3 ^ 
Mz cos Qw cos /5 — M^ cos Qw sin /? 
-/i 
-/i J 



By W„(x) we denote T'3(x) PF„2+tan6't^ -^y^ W^ni, where T^ix) is the third component 
of the weak isospin and Y{x) the hypercharge. 
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